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ABSTRACT 


It has been reported recently by McDonald ^ sd. (1976) that 1-^V 
protons may undergo considerable acceleration in co-rotating streams. 

It has also been suggested recently by Fisk et al. (19/4b> that inter- 
stellar neutral particles idiich are ionized in the solar cavity nay be 
accelerated in the solar wind, and may account for the anomalous com- 
ponent that is observed in low-energy cosmic-rays (at 'v 10 JfeV /nucleon) 
It is shown here that the particles in both of these cases could be ac- 
celerated by transit-time damping propagating fluctuations in the mag- 
nitude of the Interplanetary magnetic field (e.g. magnetosonic waves). 
The protons in co-rotatlng streams may be accelerated by transit-ti^ 
damping the small-scale variations in the field ma^itude that are ob- 
served at a low level in the inner solar system. The interstellar ions 
may be accelerated by transit-time damping large-scale field variations 
in the outer solar system. 
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The acceleration of energetic particles in the Interplanetary 
medium has become a subject of renewed Interest as the result of two 
recent observations. McDonald al . (1976), by comparing Pioneer 11 
and IMP 7 data, have found that the intensity of 1-MeV protons in co- 
rotatlng streams increases by more than an order of magnitude between 1 
and 3 AU from the sun. Although other interpretations of this d.-’ta may 
be possible, McDonald ^ al . (1976) draw thr reasonable conclusion that 
the intensity increases are compelling evidence for extensive inter- 
planetary acceleration of low energy particles. Further, the so-called 
anomalous component in low energy cosmic rays reported by McDonald et 
al. (1974) and Hovestadt aj[^. (1973) has been interpreted by Fisk, 
Kozlovoky, and Ramaty (1974b) as resulting from interstellar neutral 
particles which are swept into the solar cavity, and then ionized and 
accelerated in the solar wind. This interpretation, which is discussed 
in more detail in Fisk (1975 a, b, c) , provides a natural explanation 
for the composition of the anomalous component, which is mainly He, N, 

0, and Ne. This interpretation requires, however, that the interstellar 
ions can be accelerated in the solar wind from energies 1 keV/nucleon, 
which is the energy they obtain immediately following their ionization, 
to their ouserved energy 10 MeV/nucleon. 

A mechanism that is commonly invoked to accelerate particles in the 
interplanetary medium is a statistical Fermi process in which particles 
are resonantly scattered by Alfvcn waves with wavelengths comparable to 

the particle’s gyro-radius (e.g. Jchipii, 1971). Equivalently, the par- 
ticles are accelerated by cyclotron damping the Alfven waves. This mech- 
anism, however, has two disadvantages. First, to provide a reasonable ac- 
celeration rate for 1 MeV protons, extensive pitch-angle scattering of 
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the particles is required. For example, McDonald e]^ al^. (1976) find 
that even with unusually large Alfven speeds ( 100 km/sec), the Fermi- 

scattering mechanism can account for the observed intensity Increases 
only if the mean free path for 1-MeV protons near the orbit of earth and 

-3 

beyond is 7*10 AU. In estimating this mean free path, McDonald e^ 
al . (1976) do not include the diffusion in energy that occurs in a 
statistical acceleration process; an effect which reduces the required 
acceleration rate. Also, as is noted by McDonald al^. (1976), the 
intensity increases appear to coincide with colliding streams in the 
solar wind, where conceivably the mean free path is unusually small. 
However, the required mean free path is more than an order of magnitude 
smaller that the typical value for low energy particles of 0.1 AU 
(Ma Sung et » 1975). A second disadvantage of the Ferml-scattering 
mechanslm is that its efficiency is reduced if the Alfven waves propagate 
primarily in one direction, as is frequently the case in the solar wind. 
However, in regions where streams Interact, Alfven waves may be observed 
propagating in both directions along the magnetic field (Belcher and 
Davis, 1971). 

As is discussed in more detail in Fisk (1975 c), arguements based 
on the expected modulation of galactic cosmic rays easily demonstrate 
that the Fermi-i-cattering mechanism is not responsible for the accel- 
eration of the interstellar ions in the theory of Fisk e^ a^. (1974b). 

In this theory the particles in the anomolous component are only singly- 
charged (during their dwell time in the solar cavity, the interstellar 
neutrals can be Ionized only once (Fisk £t al. , 1974b)). Anomolous oxygen, 
then, in the energy range from 1 keV/nucleon to 10 MeV/nucleon (the range 
over which it must be accelerated) has a rigidity equal to that of, e.g., a 
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proton In the energy range from 250 keV to 1.5 GeV, l.e. equal to the 
rigidity of modulated galactic protons. Since the mean free path Is a 
function of rigidity, modulated cosmic rays and anomolous oxygen ex- 
perience the same scattering conditions. If the scattering Is sufficient 
to accelerate the anomolous oxygen by a Ferml-scatterlng process, then a 
modulation of galactic cosmic rays far In excess of allowable limits 
results . 

In this paper, transit-time damping Is examined as a possible means 
for accelerating low-energy particles in .:o-rotating streams, and inter- 
stellar ions in the theory of Fisk ^ a^. (1974b). In this process par- 
ticles are accelerated by interacting with propagating fluctuations in 
the magnitude of the interplanetary magnetic field (e.g. magnetosonlc 
waves). This mechanism overcomes the two disadvantages associated with 
the Ferml-scattering process. As will be shown, this mechanism results 
in little pitch-angle scattering. Also, unlike Alfven waves, magnetosonlc 
waves do not all have the same phase speed parallel to the mean magnetic 
field. Thus, there is an inherent dispersion to these waves, so that 
even if the waves propagate, e.g., primarily outward from the sun, a 
Fermi acceleration is still possible. 

Magnitude fluctuations in the solar wind are present only at a 

relatively low level. The dominant fluctuations result from changes in 

the direction in the field that are caused primarily by Alfven waves. 

In the inner solar system, for example, the rms change in a component 

6B^ of the field, divided by the mean field strength squared, is 
, , 2 

<( 6 B ) ^^> / B = 0.1 - 0.2. In contrast, the rms change in the 
X o 

field magnitude is i <( o B )-> / B^^ = 0.01 - 0.03 (Smith, 1974). 

TtiC m.ignitudo fJuctuatlops are presumably reduced because these variations 
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are daroped not only by energetic particles, but also, as has been shown 
by Barnes (1966, 1968a, h, I 969 ) they are Landau damped by solar wind pro- 
tons and transit-time damped by solar wind electrons. T t is hhovm here, 
however, that the observed value of yields an acceleration rate for 
protons that is sufficient to account for the observations of McDonald 
(1976). For the acceleration of interstellar ions In the outer 
solar system a larger value of will be required (n^ = 0.1). However, 
this value will be justified on the grounds that the magnitude fluctu- 
ations involved are not readily damped b‘ solar wind particles. 

Magnitude fluctuations should be a contino*’i source of energy for 
accelerating energetic particles because these fluctuations are con- 
tinuously regenerated. For example, Alfven waves decay into magnetosonic 
waves as they propagate through the spiral pattern of the interplanetary 
magnetic field (Belcher and Davis, 1971). Variations in the solar wind 
velocity result in compressions of the Interplanetary field. Further, 
a stochastic magnetic field, like the Interplanetary field, cannot be 
hydrostatically stable (Parker, 1972), and thus in continuously adjusting 
the particle and field pressure, should produce magnitude variations. 

In section 2 the acceleration rates that are required to account 
for the observations of McDonald e£ (1976), and to account for the 
acceleration of Interstellar ions in the theory of Fisk et (1974b), 
are determined. In section 3 the acceleration rates that result from 
particles transit-time damping magnitude fluctuations in the inter- 
planetary magnetic field are calculated from standard quasi-linear/ 
adiabatic theory, and are found to agree well with the required rates. 

The calculations presented in this paper provide formal justification 
for the acceleration rates used in Fisk (1975 c) . 
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2. REQUIREMENTS FOR THE ACCELERATION 


(2,1) Low-energy protons In co-rotating streams . Even with a mean 
free path as large as 0.1 AU, 1-MeV protons tend to co-move with the 
3olar wind. An estimate of the acceleration rate that is required to 
account for the observations of McDonald ejt (1976) can then be 
obtained by considering that the particles move primarily by being 
convected with the solar wind, rather than by their own diffusive motion. 

Consider that the acceleration is the result of a statistical Fermi 
proces , . or equivalently the result of the damping of a random distri- 
bution of waves. For example, the acceleration could be by the Fermi- 
scattering process (cyclotron damping), or it could be by transit -time 

damping. Such statistical acceleration can be described as a diffusion 

_ 2 

in momentum space, with a diffusion coefficient D = <(Ap) ^ /At, 

PP 

which represents the rms change in Che magnitude of particle momentum p, 
per unlc time, averaged over all particle directions. The co- 
efficient u ^quely determines the rms change in energy and the mean 
change in energy per unit time; the latter term is the one considered by 
Jokipii ( 1971 ) and McDonald e_t ( 1976 ). In a steady-state model for a 
energetic particle stream in the interplanetary medium, the particle omni- 
directional distribution function f (number of particles per unit volume of 
phase space, averaged over particle direction) behaves, then, according to the 

( 1 ) 


equation: <> o 

V 3f 13 (p^D 3f) P_ 3 (r V) 3f 

"'’a? ^ ¥ 


3r 

Here r Is heliocentric radial distance and V is the solar wind speed. 

The first term on the right side of (1) describes the effects of the ac- 
celer.atlon, and the second, the effects of adiabatic deceleration and 
convection In the solar wind (see Fisk et al., 1973). The differential 




TB 
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Intensity j , per unit Interval of kinetic energy E, Is related to f by f 


■ j/p (Forman, 1970). 
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Consider that D -Dp, where D Is a constant, independent of p 

° ° , , , -2/3 

and r. With the variable transformations r -► r , p -♦ p r , and 

then r -^ (( 7 - 23)y/3) \ (1) can be reduced to a simple 

diffusion equation of the type that is used in studying the propagation 

of energetic particles in solar flare events. Here y is the analogue of 


time in the flare oroblem, and p , the analogue of radial distance. For 

example, suppose that particles are injected into the accelerating 

region of the co-votatlng stream at a position r^ and momentum p^, such 

that r^ is much less than the position r at which observations are made, 

and p is much smaller than the observed momentum. Then, (1), written 
o 

f 

in terms of y and p , can be solved as is the flare propagation problem 
with impulsive injection of particles at the sun (Parker, 1963). The 
result in terms of r and p is 


f(r,p) “ r 


-(7-26) 


exp 


1 ( 7-2B ) 

3 2 

^ ( 2 - 6 )^ 


2 - 6 „ 

2 1 

D r 
o 


( 2 ) 


for 0 < B < 2. 

As is demonstrated below, B lies in the range 0-0.5 foi transit-time 

damping. With B 0 and with a solar wind speed of V - 400 km/sec, (2) 

- 2 -7-1 

then predicts that a value of ®pp/P * 1.5*10 sec for 1-MeV protons 

is required to account for the observations of McDonald et al . (1976). 

The intensity of these particles then increases by roughly a factor of 

ten between r - 1 and 3 AU. This value and form for D yields, from 

PP 

(2), a spectral shape at r - 3 AU of j ® E exp ( -3.5 E ], where E is in 
units of MeV. This spectral shape is not unrealistic, as can be seen by 
examining the proton intensities at different energies, reported in 
Figure 3 of McDonald e^ (1976). 
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(2.2) Interstellar Ions. Interstellar neutral particles, unlike their 


ionized counterparts, are swept into the solar cavity by the motion of 
the sun through the interstellar medium. The neutrals, which should 
consist mainly of H, He, N, 0 and Ne, are then ionized here by photo- 
ionization from solar UV and by charge-exchange with the solar wind. 

(See Fisk (1974h) for a defense of the neutral composition and 

Axford (1972) for a general review of this subject.) Once ionized, the 
inters* rellar particles, which are only singly-charged, a»*e forced to co- 
move nith the solar wind, and obtain energies in this process ''< 1 
keV/nucleon (Fisk and Goldstein, 1976). In the theory of Fisk et al. 

(197Ab) and Fisk (1975 a, c), it is assumed that a small fraction of tiie 
Interstellar Ions are then accelerated slowly but continuously in the 
solar wind, so that by the time the ions have been convected to large 
heliocentric distances (e.g. r 50 AU) , the particles have obtained 
energies 10 MeV/nucleon. Singly-charged He, N, 0, and Ne at these 
energies have high rigidities and can diffuse back into the inner solar 
system, producing the observed anomolous composition. Interstellar H, 
and any accelerated solar or galactic particles (which should be fully 
stripped of their electrons) are still at a low rigidity at these energies, 
and are excluded from the inner solar system by solar modulation (Goldstein 
e^ ^. , 1970; Gleeson and Urch, 1971). 

A rough estimate for the spectral shape of the accelerated ions in 
the outer solar system can be obtained from (1) and (2), in which par- 
ticles are assumed to move spatially by being convected by the solar 
wind, primarily, rather than by diffusing. This approximation is not 
unreasonable even at high rigidities. With relatively slow acceleration 
rates, the spatial gradients can adjust so that the loss of particles 
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from the solar cavity is no larger by diffusion than It Is by convt-ctlon 
(Fisk, 1975 c) . Detailed numerical calculations of the spectra of 
accelerated interstellar He and 0, in which the effects of diffusion, 
connection, and adiabatic deceleration have been included, can be found 
in Fisk (1975 c) . 

Anomalous He particles are observed to obtain higher energies than 

do the anomolous N, 0, or Ne (McDonald et^ al.. * 1974; Garcla-Munoz et al., 

1973). Thus, the He observations place the most severe constraints on 

acceptable acceleration theories. Reproduced in Figure 1 is the He data 

for quiet-times in 1972, as was reported by Garcla-Munoz et -il^. (1973). 

Shown also in this figure is the spectrum that is expected for modulated 

galactic He. Clearly, to produce the flat portion of the spectrum in 

the energy range from 'v 10-50 MeV/nucleon, which is the so-called anomolous 

portion, a spectrum of accelerated interstellar He at 1 AU, as is 

shown in Figure 1, is required. 

- 2 -8 —1 

With a value of D /p ■ 3*10 sec at E * 10 MeV/nucleon, and 
PP 

with B ■ 0.5 and V » 400 kra/sec, (2) yields the spectral shape at r ■ 50 

AU that is shown in the upper curve in Figure 1. The normalization of 

this curve is arbitrary since it depends critically on the exact number 

of Interstellar ions that are actually accelerated. This number is 

-4 

expected to be only a small fraction ( ''^ 10 ) of the number that are 

ionized in the solar cavity (Fisk, 1975 a, c). The modulation of the 

interstellar He, which is not excessive since the particles are singly- 

charged, should harden the spectrum as the particles propagate from r ■ 

50 AU back Into the Inner solar system. Thus, this value of 5 should 

PP 

be able to produce the required accelerated spectrum at 1 AU. 

Note that the acceleration rate required for the Interstellar ions 
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is roughly equal to the rate required for protons In co-rotatlng streams. 

-2-15 - 2 

It was assumed that D /b « p * . Thus, at 1 MeV/nucleon D /p *■ 1.7* 

PP PP 

lO'^sec”^ for Interstellar He, compared with D^p/p^ * 1.5*10 ^sec ^ for 
the co-rotating protons. This comparison, however, is not necessarily 
meaningful since the acceleration in these two cases can occur under 
substantially different conditions. The acceleration in co-rotating 
streams is a local phenomenon in the inner solar system, which appears 
to occur near regions where solar wind streams interact. Interstellar 
ions are accelerated over much larger dimensions. To reach the required 
high energies of 10-50 MeV/nucleon, interstellar He must be accelerated 
continuously as it is convected out to distances of r 50 AU. 

3. ACCELERATION RATES 

^ 'y 

(3.1) The formal calculation. The acceleration rate D /d" thnt resultK 

PP 

from particles transit-time damping fluctuations in the magnitude of 
the Interplanetary magnetic field can be calculated by applying the 
standard procedures of quasi-linear /adiabatic theory, as has been done 
in similar contexts by Barnes (1968a). Consider that the magnitude fluc- 
tuations result from the superposition of fast-mode magnetosonic waves. 
This mode is more likely to be present in the interplanetary medium than 
is, e.g., slow-mode waves (cf. Hollweg, 1975). In the hydro-magnetic 
limit the fluctuating magnetic field of a fast-mode wave must lie in a 
direction that is both normal to the wave vector k and normal to a unit 

vector 8, which lies in the B x k direction; B is the mean magnetic 
1 -o — -o 

field (e.g. Lee and Vtflk, 1975). The fluctuating magnetic field at 
position £ and time t is then 

B'(r,t) • A()j) (K ^ 

k 
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Here, the function A fk) is defined such that A (k)A* Qc) = C (k) , 
where the star denotes complex conjugate and 

rd\ C (k) - <(6B)^> (4) 

•00 ““ 

2 

with <(6B) > the rms change in the field magnitude, associated with the 
waves. From Maxwell's equations, the fluctuating electric field Is 

E'(r,t) - /*d^k/*du) exp li(k*r - wt)! ojA(k)fi, (5) 

•“ " 00 •OD — • ““ ^ 

ck 


where c is the speed of light. In this calculation only normal modes 
of the plr.sT^ \re considered, which requires that the frequency w satisfies 
the relationship 

J - k^ U/2(V^ + vh > 1/2 I (vj V V - 4V^ k?,/k^]^^^} (6) 

AS AS AS 

Here V. and V_ ary the Alfven and sound speeds, respective!' , In the 
A S 

solar wind; k„ is the component of k parallel to B^. The turbulence is 
also taken here to be homogeneous. 

The fluctuating electric field will change the magnitude of the mo- 
mentum of energetic particles. By following standard quasl-1 inear/ad 1- 
abatic theory (cf. Goldstein , 1975), this change can be described 

by a diffusion coefficient. 




2 


PP 


q / dt / d k exp Il(k*f - u t)] C(k)y 

« — “O — opiip 


2 0 


Here q is particle charge and is a unit vector in the direction of 
particle momentum £. The tilde denotes streamed variables, which are to 
be evaluated along the trajectory of a particle spiralling along 
Streamed variables are then functions of t. 

After considerable manipulation, (7) can be reduced to 
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D 

_E£ 

2 

P 


_l /"dk„ /"k dk, C (k,. k„. ♦) J u'^ I 

2 2 0 - 0 - - j 

TO k 


( 8 ) 


00 2 
where I * ti E 6 (k„vy-a) + mil) J ' (u'k r ) 
m— « o m i g 


(9) 


Here k is defined in a cylindrical co-ordinant system, with the axis of 
symmetry k„ and k, are measured parallel and normal to respec- 
tively, and is the azimuthal angle. Particle gyrofrequency in the 

mean field B is and particle Larmour radius is r . Particle speed 
-o g 

2 1/2 

is V and pitch angle y, with y' = ( 1 - y ) . The function J'^ (z) 

is the derivative with respect to z of a Bessel function of the -first 
kind and order m; 6(z) is the delta function. 

The diffusion coefficient for pitch angle scattering by magnetosonic 
turbulence, (•*<(Ay) ^>/At ) , can be found by the same procedure to be 

2 

% = /“dk„ /“kjdk^ c (k,, k„, I (10) 


The scattering rate is roughly the inverse of the mean collision 
time, or equivalently it yields a mean free path 


A reasonable form for C (k,,k„, i|;) is 
2 2 

^ 2n B r(v + 2) X X X 

C (k,, k„, ij;) ^ o X y z 

r(v-l/2) 


2 2 2 2 
(X‘ kf, + X, kp 

Z JL 1 

j. 2 . ,2 , 2.V+2 
(1 + X k„ + X k ) 

Z JL J. 


( 11 ) 


where X = X ^ cos^ 4; + X ^ sin^ 4< and = <(6B)^>/B Here, X , 

1 X y o z 

\ and X are the correlation lengths for the turbulence in the 
y X 

direction parallel to ^ and the two orthogonal directions normal to 
respectively; I‘(z) is the gamma function. The index v must exceed 0.5 
so that the energy density in the field fluctuations is finite. If the 


fieln is sampled in, e.g., the direction along this form for C yields a 
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one dimensional power spectrum for fluctuations parallel to ^ which ap- 
proaches a constant as k„--0 and declines as k„ for large k„. With v 
in the range 0.5 - 1.5 this power spectrum resembles observed power 
spectra (*Toklpii and Coleman, 1968; Sari and Ness, 1969, 1970). 

The behavior of C near k = 0 given in (11) is chosen so that the 
derivative of the correlation function in Ic-space are well -defined at 
this point. Equivalently, the moments of the correlation function in r- 
space are required to be finite. Dependences near k = 0 in which C de- 
creases less rapidly with and for which the derivatives are not well 
defined, are also possible. In practice, the behavior of C near k = 0 
is relatively unimportant since it influences the acceleration rate of 

only quite high energy particles (i.e. r >> X , X , A ). The behavior 

g X y 2 

of C near k = 0 in (11) is in fact a conservative choice since it yields 
less acceleration for high energy particles than do other acceptable de- 
pendences . 

It should be noted that in deriving (8) and (10) the particle dis- 
tribution function was assumed to be gyrotropic about This as- 

sumption is technically incorrect for turbulence that is not axial-sym- 
metric about B , i.e. it is incorrect when C varies with However, 

— o 

deviations from gyrotropy should occur only on times long compared with 
the time required for the integral over t in (7) to approximate closely 
its limiting value. In practice, the gyrotropic approximation should 
introduce no significant errors, even under the extreme circumstances 
considered below. 

(3.2) An illustration example. Consider the case where v > 1 and r 

g 

<< X^, X^, and X^. This example is instructive in that it readily 
illustrates the physical processes involved in transit-time damping. 
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To a high accuracy only the m = 0 resonance need be considered in 
2 

calculating D /p from (8). With v > 1 and r << A A , the Bessel 
PP g X y’ 

function in this resonance can be well-approximated in ( 8 ) by J '(ii'k r ) 

o i g 

^ -p'k,r^/2. Further, to the accuracy required here it is reasonable to 
use the approximation that * + u, where u is a constant phase 

speed. The weak dependence on tp in the dispersion relation for fast- 
mode magnetosonic waves given in (6) is thus Ignored. With this last 
approxlnation the resonance condition prescribed by the delta function 
in (8) requires (for m = 0) that 


k,,viv = w s: + uk 
o — 


( 12 ) 


which for energetic particles (p v >> u) yields 


+ u k^ 
pv 


(13) 


or thus k„ << kj . In (8), men, terms containing k„ can be neglected 

relative to terms in k^, provided that the ratio of A^ to A^ and A^ Is 

such that A u/(yv A ) and A u/(uv A ) are also small. 

2 X z y 

With these approximations (11), (9), and (8) give 


D 2 ,4 2 

V 

„2 ' vlu 


■ ;(v+2) AAA 
X y z 


r (v - 1/2) 


2 


/ "dtJ^A, / dk.k^ (1 + A^ k^)“'^"^ , v>l 
0 ■ 0 " * ■ 


(14) 


which can be Integrated to yield 


D 2 ,4 2 
PP u u' n 


r(v-l) 

6 r(v-i72)" 


X 2 A 2 

X y 


( 15 ) 


The terms in (15) can be readily identified, by the following 
argument . 
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A particle that Interacts with a propagating fluctuation In the 
magnitude of the magnetic field experiences a change In the magnitude of 
Its momentum, which can be expressed as 


P 


u„Av„ 


(16) 


Here u„ Is the phase speed of the fluctuation along the mean field, and 

Av„ Is the change in v parallel to the mean field, as viewed from the 

frame moving with u„. In deriving (16) use is made of the fact that the 

energy is constant in the frame moving with u„. 

As can be seen from (14) and (13), with v>l the fluctuations that 
2 

contribute most to D /p are those for which k„ = l/> ' << 1/' and k 

pp " z z ' 

2 2 1/2 

= (1/A + 1/A ) . For such fluctuations u,, has a magnitude u„ 

X y 

2 2 1/2 

o- uA ' ( 1/A + 1/A ) . Since there is no electric field in the 

7 X y 

frame moving with u„, Av„ results from 8B/3z, the gradient in the 

magnitude of the field, along B^, which is produced by the fluctuation. 

With r << A , A and A 
g X* y 2 


dv„ 

dt“ 


,2 2 

VI V 


3Z 


(17) 


or 


A I .22^ 
Av„l ^ yi V nAt 


(18) 


Hence 


P 


2 

^ vj* unAt 


2 2 

A A 
X y 


1/2 


(19) 


Here, At ''' A^/(|u|v) is the time over which the particle interacts with 


a single, coherent fluctuation. 
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Clearly, to within a numerical constant D /p in (15) is 'V' 

PP 


2 

(Ap/p) /At where Ap/p is given in (19). 

In transit-time damping, then, particles are accelerated by inter- 
acting with a series of randomly moving gradients in the magnetic field. 
Although Av„ is considered to be small in transit-time damping, this 
process is closely analogous with a classic Fermic process in which par- 
ticles are accelerated by interacting with moving magnetic mirrors. 

2 

It should be noted that D /p in (15) varies inversely as the 

PP 

square of the correlation lengths normal to ^ and directly as the cor- 
relation length along Thus, the acceleration rate tends to be large 

when the characteristic length-scales of the turbulence are relatively 
small in general, or when the turbulence is elongated such that >> 

X and/or X . 

X y 

It should also be noted that particles interact here with long- 

wavelength waves. Thus, the magnetic moment is a constant of the motion, 

as was assumed in deriving (17) and (18). Transit-time damping, then, 

can Increase the energy of a particle only in motion parallel to or 

equivalently transit-time damping decreases the average particle pitch- 

’4 

angle. The acceleration rate, however, varies as \i . For the accel- 
eration to remain efficient, then, pitch-angle scattering must iso- 

2 

tropize the particle distribution on a time-scale of at least p /^pp* 

The pitch-angle scattering that normally occurs in the Interplanetary 

medium (X - “v O.IAU) is generally more than sufficient for this purpose, 
mrp 

(3.3) A more raalistic example . One-dimensional power spectra that are 

observed in interplanetary space generally fall off less steeply than 
_2 

k for large k, which implies that v<l (Jokipli and Coleman, 1968; 

Sari and Ness, 1969, 1970). With such values for v, the integral over 
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k, in (14) does not exist. The approximation J '(y'k, r ) - -p'k, r /2, 
i o i R 1 R 

which was used in deriving (14), is not valid at large values of k^, 

and with v<l causes the Integral to diverge. 

2 

An accurate value for ^pp/P when v<l can be obtained by integrating 


(8) and (9) numerically, as is dene below. An analytic approximation to 

Dpp/p can also be obtained by noting that the Bessel functions in (9), 

rather than causing (8) to diverge, actually truncate the integral over 

k, at values of k, > l/u*r . This truncation is to be expected because 
i if? 

the effects of fluctuations that have scale-s.^zes small compared with 


p'r^ should average to zero, ever a gyro-perlod. To a crude approxl- 
2 

nation, then, D /p can be estimated by letting J *(u'k,r ) = - p'k r /2 
pp ^ ^o'-lg ±R 

and integrating only over the range 0 < k^ < 1/p'r^, or 


D 2 '4 2 

pp . u p n 


r(v+2)x >. X 

X y z 

6iT^^^r(v-l/2) 




.1/p’r 


2.-V-2 


R dk,k^ (1+X^ k^) ^ 


With p'r << X , (20) yields a substantially larger acceleration 
R i 

rate than does (14). This Increased acceleration can be attributed to 
the smaller length-scales associated with the turbulence when v*"!. As 
can be seen in (19), for example, the momentum change of a particle inter- 
acting with a wave varies directly as the amplitude of the wave and in- 


versly as the wavelength normal to With v>l, this ratio is simply 

2 2 - 1/2 

6B over the perpendicular correlation length (1/X + l/X ) . With 

X y 

v<l, however, there is sufficient energy in the smaller wavelength fluc- 

2 2 1/2 

tuatlons to yield a ratio large compared with 6B(1/X^ + 1/^^ ) 

2 

It should be noted that D /p in (20) is now a function of r . 

PP g 

For r '<'< X and X , it can be shown that D /p^ varies as r 
g X y pp R 


( 20 ) 
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(v<l)« The momentum dependence Is “ p 

2 

When u'r >>\ and X , D /p in (20) can be integrated to yiidd 

R X y PP 


2 '4 2 
u y n 

1 /2 

IT '"^r(v+2) 

X 



v|y I 

36 r(v-l/ 2 ) 

z 

,2 

.2 


_ 


LjS JLl 



( 21 ) 


Upon comparing (21) with (15), it can be seen that the acceleration is 

inefficient when y’r >>X and X . 

R X y 

(3.4) Comparison with required acceleration rates . To compare the 

acceleration rate for transit-time damping, with the required rates 

discussed in Sections 2.1 and 2.2, (8) and (9) must be averaged over 

_ 2 

particle pitch-angle to form ®pp^P • *»traight forward average 


5 - f]du D 

PP -1 PP 


( 22 ) 


however, does not seem to be a meaningful procedure. As can be seen in 

2 

(15), (20) and (21) ^pp/P becomes large as y->0, and thus in (22) 

depends strongly on the behavior of D at small values of p. Unfor- 

PP 

tunately, at small y quasi-linear/adiabatic theory is unlikely to be 
correct (cf. Fisk e;t ad. , 1974a). 

It should be noted, however, that the Integral in (22) does exist, 

despite the apparent divergence of as y->0. The approximations 

(yvX )'<^'l and X u/(yvX )<<1, which were used in deriving (15), (20) and 
X z y 

(21) fall for small y, and in fact if computed accurately, » 0 at 

y - 0. 

Rather than performing an average, the expedient and conservative 

procedure that is followed here is simply to take D = D , where D 

PP PP PP 

is evaluated at a characteristic value of y > 0.3. 

(3.4.1) Low-energy protons in co-rotating streams . The co-rotating proton 
streams reported by McDonald e^ a^. (1976) appear to coincide with stream- 
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stream interaction regions in the solar wind. In such regions it is 

2 10 
conservative to take n ■ 0.01, u ■ 75km/sec and X "X * 2*10 cm. 

X z 

2 

This value for n is typical of interplanetary conditions iu the inner 
solar system (Smith, 1974). Similarly, in specifying u it is assumed 
simply that the Alfven speed and ion thermal speed are equal at about 50 
km/sec, which is again characteristic of average conditions. The 
correlation lengths used here are equal to the value found by Fisk and 
Sari (1973) for small-scale propagating waves in the solar wind. This 
correlation length is smaller than the value of 2*10^^ cm given by 
Jokipii and Coleman (1968); however, Fisk and Sari (1973) argue that the 
length 2"10^^cm applies only to large-scale tangential discontinuities 
in the interplanetary field. Further, it seems reasonable to use a 
small correlation length for an interaction region, in which the turbulence 
may be generated locally. 

Upon assuming that the particles are protons, that * 5*10 ^Gauss, 

and that v = C.75 and p » 0.3, (8) and (9) have been integrated numeri- 

2 

cally with the above values for n u, and X , X and X . Several 
^ X y z 

resonances were computed, but again only the m “ 0 resonance is Impor- 

2 

tant. Shown in Figure 2 is ®pp/P ^ function of energy. It can be 
- 2 -7 

seen here that D /p*" = 1.2*10 at 1 MeV, which is close to the ac- 
PP 

celeratlon rate that is required to account for the intensity increases 


reported by McDonald (1976) (cf. Section 2.1). 

- 2 -2 

It should be noted in Figure 2 that ^pp/P varies roughly as p 

- 2 -1.5 

This dependence is intermediate between D /p ® p * , which holds when 

PP 

- 2 -7 

r^<<X^,X^, and ^pp/p “ P which, from (21), is the case when • 

_ 2 —2 6 

The dependence D /p “ p yields B=0(D -Dp), which is the 

PP PP * 

value used in the discussion in Section 2.1. 
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Also plotted In Figure 2 Is the scattering mean free path 

'V' v/D . where D Is calculated from (10) by using the same parameters 
PU WP 

2 

as were used in computing ^pp^P * most important contribution in 

this case comes from the m ■ 1 resonance, l.e. the scattering is due 
mainly to fluctuations with scale-sizes comparable to the particle gyro- 
radii. The mean free path here is larger than the typical value of 0.1 
All, and thus acceleration by transit-time damping does not require or 
imply extensive pitch-angle scattering. Presumably, pitch-angle scat- 
tering in the interplanetary medium is the result primarily of particles 
interacting with Alfven waves, which are not considered here. 

(3.4.2) Interstellar ions . As was discussed in Section 2.2, the accel- 
eration of interstellar He ions in the solar wind is assumed to occur 
continuously so that as the particles are convected to large helio- 
centric distances they obtain energies of ^ 10 MeV/nucleon. In the 
outer solar system, the magnetic field is relatively weak and thus the 
gyro-radius of a 10-MeV/nucleon particle is large. For example, if the 
field strength falls off as 1/r, which is the case if the field pattern 
is that of an Archimedes spiral, then at r * 50 AU, 'v- 10 ^ Gauss, and 
r^ for a singly charged, interstellar He ion at 10-MeV/nucleon is '^-O.IAU. 
As can be seen from (20), however, acceleration by transit-time damping 

is inefficient when r >>X ,X . Thus, acceleration in the outer solar 

g x’ y 

system can result only from the interaction of particles with large- 
scale field structures. 

In addition to small-scale turbulence, such as is used above in 
accelerating co-rotating protons, stream-stream interactions in the solar 
wind may also produce the required large-scale variations in the inter- 
planetary magnetic field. As a fast stream overtakes a slow one, the 
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. . I 


solar wind density, pressure, and magnetic field tend to be compressed 
throughout the interaction region. As is discussed in detail by, e.g.. 
Gosling ^ al. (1975), within the first few AU from the Sun these com- 
pressions become enhanced with increasing heliocentric distance, and 
eventually steepen into shocks. The compressions at this point occupy 
isolated regions of space, a result which is compatible with the Pioneer 
10 observations (Smith and Wolf, 1975). The shocks, however, propagate 
relative to the solar wind and eventually should cause the compressions 
to spread and thus to disipatt . Tliis spreading could conceivably cause 
adjacent compression regions to merge with each other. Thus, at large 
heliocentric distances, e.g. r>10 AU, the compression regions could 
conceivably dislpate into a homogeneous distribution of propagating 
large-scale fluctuations in the magntiude of the interplanetary field. 

These large-scale fluctuations should have a highly elongated shape 
in the outer solar system; a configuration, which as was discussed in 
Section 1.2, should make them efficient accelerators. The solar wind 
expands in the azimuthal and polar directions about the Sun, but does 
not expand appreciably in the radial direction. With the field in an 


Archimedes spiral pattern, X is measured in essentially the azimuthal 

z 

direction at large r; is taken to lie in the polar direction, and 


then X in the radial direction. Thus, X and X both increase in pro- 
y z X 

portion to r, with X^ held constant. 

2 

Consider that n ■ 0.1, u ■ 73km/sec, v ■ 0.75 and u ■ 0.3. The 


correlation lengths are taken to have the above radial dependences, and 


at r ■ 1 AU to have the values X ■ 0.2 AU and X ■ X ■ 0.4 AU. The 

y X z 

mean field is assumed to be B ■ 5*10 ^ Gauss at r ■ 1 AU, and to de- 

o 

crease as 1/r. With these values and upon assuming that the particles 
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are singly charged He ions, ^pp^P been calculated by Integrating 

(8) and (9) numerically. 

2 

Plotted in Figure 3 is ®pp/p ^ function of energy at various 

heliocentric distances. It can be seen here that at r ~ 50 AU and at an 

— 2 -8 ~1 

energy of 10 MeV/nucleon D /p = 3*10 sec . As was discussed in 

PP 

section 2.2, this rate is sufficient to account for the acceleration of 

interstellar ions in the theory of Fis/ £t (1974o) and Fisk (1975c). 

At energies above 1 MeV/nucleon, varies roughly as p 

Thus. 5 ■ U which is the functional form used in Section 2.2. 

PP 

2 . 2 
The value of n =0.1 used here is larger than the value of n * 

0.01 - 0.03 which is observed in the inner solar system (Smith, 1974). 

2 

Further, n is taken to be constant here, independent of r, and thus no 

2 

consideration has been given to the possibility that ri" may decrease 

with r as the waves do work on the solar wind (Hollweg, 1974). The 
2 

constant value of h * 0.1, however, could be justified in the outer 

solar system, since large-scale elongated fluctuations may not be readily 

damped by solar wind particles. As can be seen in Figure 3, at large 

heliocentric distances the acceleration rate for 1-kev/nucleon ions 

becomes small. In this region, the fluctuations have become sufficiently 

elongated so that for these low-energy ions the ratio X u/(X |u|v)^>l. 

z y 

Thus, the transit-time of a particle over a correlated fluctuation 
(X^ /IpIv) Is long compared with the characteristic time for the fluc- 
tuation to change ^^/u. The particle gains only a small portion of the 
energy that the gradient in the field strength can impart, before the 
gradient changes, and thus the acceleration is reduced. The speed 
of, e.g., solar wind electrons is comparable to the speed of 1-kev/ 
nucleon ions, and thus the acceleration of, or the damping by, solar 
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I 


I 


wind particles should be slmilaily reduced. 

The Interplanetary magnetic field should lie In the azimuthal di- 
rection in the outer solar system only in some average sense. Sig- 
nificant local deviations from this direction are expected. In cal- 
culations of the behavior of the scattering mean free path with helio- 
centric distance, these local deviations have an appreciable effect. 

For example, Vblk a^l. (1974) found that if the field lies strictly in 
the azimuthal direction, the spatial diffusion coefficient increases 
rapidly beyond earth and becomes Imcompatlble with modulation obser- 
vations. In a subsequent paper, however, Morfill et al^. (1975) found by 
including local deviations that the diffusion coefficient is consistent 
with, for example, measurements of the cosmic-ray gradient from Pioneer 
10 . 


In contrast, local deviations in the field direction do not appear 
to have an appreciable effect on the acceleration rate from transit- 
time damping. As can be seen from (19) the momentum change of a par- 
ticle interacting with a propagating fluctuation In the field magnitude 

depends on the nvignltude of the k-vector of the fluctuation, not on its 

2 

direction relative to B . In the case considered in (19), Ikl (l/> 

2 1/2 

+ ) • Thus, the field can undergo small-scale changes in di- 

rection within a correlated segment of turbulence (these directional 


changes could be the result of, e.g., 
affecting the total momentum change, 
the local field direction can also be 


Alfven waves) without appreciably 

2 

The insensitivity of D /p to 

PP 

seen in (14) and (20). The in- 


tegral in these expressions depend only weakly on k„. 
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CONCLUDING REMARKS 

It was shown in this paper that acceleration by transit-time 
damping can account for the intensity increases In co-rotating proton 
streams that have been reported by McDonald et al. (1976). The small- 
scale fluctuations in the magnitude of the Interplanetary magnetic field 
that are observed in average conditions in the inner solar system appear 
to yield a sufficient acceleration rat.’. It should be noted here, 
however, that in these calculations the waves were assumed to be fast- 
mode magnetoeonlc waves and the correlation length was taken to be 
reasonable”, but small. At times when the fluctuations are cauF.ed by 
other modes or the correlation length is larger, the acceleration rate 
could be less. Indeed the coincidence of the events observed by McDonald 
£t . (1976) with stream-stream interaction regions in the solar wind 
may result because fast-mode magnetosonic waves and small correlation 
lengths are characteristic of such interaction regions. 

It was also shown here that interstellar particles which are ionized 
in the solar cavity can be accelerated to energies > 10 MeV/nucleon by 
transit-time damping large-scale variations in the magnitude of the 
Interj etary mangetic field in the outer solar system. Ttiese inter- 
stellar particles can thus account for the anomolous component observed 
in low-energy cosmic-rays, as has been proposed by Fisk e^ (1974b). 

It is .suggested here that the required large-scale field variations 
could be the remnants of stream-stream interaction region in the solar 
wind. Of course, this is only speculation, and only with direct ob- 
servation or with more detailed models of the evolution of solar wind 
streams will it be possible to determine whether the required fluc- 
tuations exist. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


The differential intensity of accelerated interstellar 
He that is required at 1 AU. This intensity when com- 
bined with the expected flux of galactic He produces 
the anomolous or flat portion of the observed spectrum. 
Also shown is the intensity of accelerated interstellar 
He at 50 AU, which is predicted from (2). 

_ 2 

A plot vs. kinetic energv of the acceleration rate D /p 

pp 

that results from protons transit-time damping fluc- 
tuations in the magnitude of the interplanetary magnetic 
field. The parameters used in calculating this curve 
should describe the small-scale fluctuations in the field 
magnitude which occur in the inner solar system, in 
stream-stream interaction regions of the solar wind. 

Also shown is the characteristic mean free oath for 

pitch-angle scattering which results from the magnitude 

_ 2 

fluctuations that are used in calculating D /p . 

PP 

A plot vs. kinetic energy, at various radial distances, 

2 

of the acceleration rate D /p that results from 

PP 

singly-charged interstellar helium ions transit-time 
damping fluctuations in the magnitude of the inter- 
planetary magnetic field. The parameters used in 
calculating these curves may describe large-scale varia- 
tions in the field magnitude in the outer solar system. 
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